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Abstract 

The theoretical investigation of the cluster de Haas - van Alphen (dHvA) oscilla- 
tions in three-dimensional systems performed for the first time. Applying a three- 
dimensional oscillator model to systems with electron numbers 10 < N < 10 5 we 
predict distinctive size effects: the dHvA oscillations can be observed only within a 
certain temperature range determined by ./V; the lower size limit for 7V is ~ 20; the 
amount of the dHvA oscillations is reduced with decreasing N which is accompanied 
by stretching the period of the oscillations. 
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1 Introduction 



With increasing a uniform magnetic field applied to a macroscopic body its 
electron characteristics can reveal low-temperature oscillations manifesting 
themselves e.g. in the de Haas - van Alphen (dHvA) or Schubnikov-de Haas ef- 
fects [1]. Oscillations in magnetization of metal thin films, the period of which 
increases with decreasing thickness, were theoretically predicted in Refs. [2,3]. 
The Shubnikov-de Haas oscillations were experimentally observed in antimony 
plates (thickness is of order of 10 2 -j- 10 3 nm) [4]. Now the dHvA effect is 
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rather intensively studied in two-dimensional mesoscopic systems (e.g. [5,6,7] 
and works cited therein). 

However the possibility of the dHvA oscillations in atomic clusters has not 
been established so far. In this Letter we predict the properties of the dHvA 
oscillations in three-dimensional (3-D)-mesoscopic systems in the absence of 
the disorder and interactions assuming that noninteracting electrons are con- 
fined in an oscillator well. We find that in 3-D-mesoscopic systems, where 
electrons confined in all three directions, along with the dHvA oscillations 
there exist much more types oscillations than in bulk metals, thin films and 
two-dimensional mesoscopic systems. We consider the finite size effects in the 
dHvA oscillations in clusters with electron numbers 10 < iV < 10 5 in a wide 
range of arbitrarily oriented magnetic fields and find that number of the dHvA 
oscillations is reduced with decreasing N and these oscillations appear at a 
temperature which depends on N. From our calculations follows that grains 
of materials with small the Fermi energy and the effective electron mass are 
best suited to observe the cluster dHvA oscillations. 



2 Model 



The harmonic oscillator model is widely exploited in describing structure of 
quantal systems including mesoscopic ones [3,8,9,10]. This is caused by the 
possibility of obtaining simple analytical solutions for oscillator energy eigen- 
values and wave functions. Electron motion in the anysotropic 3D - oscillator 

— * 

potential and a uniform magnetic field B is described by the Hamiltonian 
H = _E i^(Pu ~ -Kf + ^Irl + f^xWU (1) 



(A = rxB /2, m* is the effective electron mass, Q u are oscillator frequencies 
at B = 0, g is the effective Lande factor, fi* B = eh/2m*c, a v is the Pauli ma- 
trix). Eq. (1) being expressed through oscillator quanta, bj^, b v (//, v = x,y,z) 
is a bilinear form of these operators (b'b, b'b\ bb). We use the boson Bo- 
golubov transformation to remove terms with two creation and annihilation 
bosons and at the same time to diagonalize the rest part, i.e. to reduce it to 
^Wab^ba, a = 1,2,3 where W a are the eigen frequencies of new oscillator 
quanta in arbitrarily oriented field. As known both fermion and boson Bogol- 
ubov transformations result in equations quadratic in eigenvalues. Therefore 
we arrive at the cubic equation in W 2 since Eq. (1) deals with three spatial 
degrees of freedom. Relative simplicity of the Hamiltonian (1) is embodied in 
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a rather compact equation for eigenvalues: 

II (nl-W')-W 2 £ uj 2 M~W 2 )=0, (2) 

v=x,y,z u=x,y,z 

— * 

uj is the cyclotron frequency, uj = 2fi* B B. 

Three solutions of Eq. (2), W + , W_, Wq, determine the single electron energies 
(hereafter h — 1) 

e(ra+, n_, n , uj) = W + {n + + 1) + Vy_(n_ + 1) + W (n + -) ± |w. (3) 

n±, no are integers. The spin contributes to the last term in Eq. (3). Solutions 
W±, Wo vary with increasing uj: W + grows, falls down and Wq depends 
mainly on the direction of B. If the field is directed along one of the symmetry 
axes of the spheroid, e.g. uj z — cu, u> x — cu y — 0, W = Vt z and W± are found 
straightforwardly. For axial symmetry (Q x = Q y ) the solutions coincide with 
those given by Fock [11]. In high fields (uj 2 ^> J2u solutions of Eq. (2) gain 
the values: 



W 2 + 



UJ 



D 

~~ 2' 



Q 2 



uJ 2 u n 2 JuJ 2 ; D 



w 2 



UJ 



w 2 



n 2 - 



(4) 



To model single-electron levels we adopt the relation [12] between semiaxes 
of ellipsoidal 3D - clusters (a x , a y , a z ) and their oscillator frequencies 



£f 



^X^X QyQ>y Vt z (X z ^ f2 /O ATM /S ' (^/ 



a x a y a z = R 3 , R = r s iV 1//3 , r s is the radius of a sphere with the volume V/N. 
Q 3 = Q x Q y Q z = W+W-Wq is an invariant. ef is the Fermi energy. Eq. (5) is 
true for comparable semiaxes and small frequencies, £l x ,y,z ^ £ f- 



3 Results of calculations 



Our calculations of the electronic susceptibility (%) and heat capacity (C) 
are performed for ellipsoidal clusters. Values of C are presented here in order 
to demonstrate the temperature averaged level density near the Fermi level, 
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which is practically proportional to C. Calculations involve so wide single- 
electron spectra that their enlarging has no effect on results presented below. 



At very low temperatures (T < A, A is a mean level spacing) x an d C are 
strongly oscillating functions of Sp/u revealing no evident periodicity v.s. the 
field (Fig. 1). Numerous peaks in \ an d C result from crossings of the Fermi 
level and an upper level nearest to it that gives at low temperatures a maxi- 
mum in x an d simultaneously a minimum in C [13] as shown in the insert in 
Fig. 1. These low temperature oscillations, being absent in bulk metals and 
thin films, are caused by the discreteness of the electronic level spectrum. In- 
creasing temperature drastically changes this picture damping high frequency 
oscillations. Starting with a temperature T start one can observe the emerging 
of regularities near the quantum limit, ep/oo ~ 1 (Figs. 2). The final sup- 
pression of high frequency oscillations occurs at an optimum temperature T opt 
(Figs. 3-6) exposing a maximum amount (n max ) of periodic oscillations. They 
are analogous with the dHvA oscillations in bulk metals and are determined 
by the common physical reason: the crossings of the Landau levels with the 
Fermi surface in a strong magnetic field. Further increasing the temperature 
(T 3> T opt ) leads to damping of the cluster dHvA oscillations (Fig. 3). 

To interpret these results we have found analytical expressions for the 3D 
oscillator level density po{s) an d twice integrated level density P2{z) which 
determine the grand canonical C and x at a fixed particle number N 



C = k, 



X 

XL 



i(3 2 J dep (e) 



86! 



de 



d 2 p 2 ( e) 
du 2 



- Po{e) 



'ay 



$(£); 



$(e); 



(6) 
(7) 



1 + C 



, /3 = l/h B T, 



Xl = —p 2 3N/2VeF, V being the volume of the system. A is the chemical 
potential practically equal to Bp and weakly oscillating with u in a wide region 
of ui up to uj ~ e F where A markedly decreases. 



Pn\e) 



f 

J Q n 



167r(i) 



3+n 



W + q W_q W q 
sin sin sin 



(8) 



The integration path (£) envelops all poles including q = 0. The residue of 
p n (e) in q = gives (after integrating over e in Eqs. (6), (7)) smooth functions 
of u for C and x {C s , Xs) 
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C s = k 2 B ir 2 TN 



Xs/ | XL 



12a=x,y,z ^ 



12A 2 
3 



- 2 + -^\T F , 



A 



(9) 
(10) 



These expressions are true at T <C ep- The first negative term in xs/ I Xl \ 
is the "3D-oscillator" Landau diamagnetism and the second (~ g 2 ) term is 
the "3D-oscillator" Pauli paramagnetism. Their absolute values are two times 
more than those for free electrons, nevertheless their ratio at g — 2 is equal to 
the well known value (xl/xp — ~ 1/3) ■ Cs and xs are those middle levels on 
which all types of oscillations are superimposed. At high enough temperatures 
(but T/ef <C 1) the most part of oscillations are suppressed and C and x 
become equal to Cs and xs- Fig- 3 shows that at 4T opt the oscillations in the 
system with N = 10 5 disappear in weak fields (ef/oj > 25) and x becomes 
equal to xs- 

Oscillating components of C and x (which are determined by residues in points 
q 7^ 0) can be divided into several constituent parts. The large scale oscillation, 
which are denote as the cluster dHvA oscillations by analogy with bulk metal 
oscillations, are described by functions cos(27rAn + /W + ) cos(7r guon + /2W + )n+ 2 . 
For the dHvA oscillations frequencies W+, Wq must not be multiple, i.e. 
W+/W- ^ n + /n^ and W + /W ^ n + /n . Since A < e f and W+ ~ uo for 
uj ^> Q the period (t+) of the main "tone" of these oscillations (n+ = 1) on 
the scale e f /oo is almost equal to 1: 

cos^ttA/WM = cos ( 2tt— — ) , t+ = £ JL^±, 

\ t+ UJ J A UJ 



though with decreasing ui the period t+ is stretched, and the smaller is the 
particle number the stronger is the stretching. In Fig. 4 this stretching is quite 



evident for iV = 10 . For instance in a simple case ui z — uj, ui x 



cu,, 







W+/uo ~ 0.5 



-2/3 



+ 1 



1/2 



+ 1 



The role of the effective electron Lande factor, entering into cos(iigujn + /2W + ) , 
turns out to be similar to that in the bulk dHvA oscillations [1] (Fig. 5): At 
g = 2 the oscillations are paramagnetic while at g — they are diamagnetic. 
At g — 1 the main tone (n + = 1) disappears in strong fields that leads to 
lowering amplitudes and doubling frequencies. Thus experimental studying of 
the cluster dHvA oscillations could give information about cluster g-factors. 
Functions in Fig. 1-4 are calculated with g = 2. 
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The temperature damping factor D x (T,n a ) of the cluster dHvA oscillations 
(n a = n + ) and the other types of oscillations considered below has for the 
susceptibility the same form as for the bulk metal oscillations [1] 

D x (T,n a ) = -^; x a = 2n 2 k B Tn a /W a , (11) 



D X (T, n+) -fact or causes attenuation of the dHvA amplitudes with decreasing 
uj as seen e.g. in Fig. 3. At T <C uj the temperature damping of the dHvA 
oscillations is inessential that allows summing over n+ to be performed. It 
leads to parabolic segments in x an d C which are well illustrated in Fig. 2 by 
the example of x- 

The temperature damping of the heat capacity oscillations is governed by 
D c (T,n a ) 



D c (T,n a ) = 



3x a 
sinh x a 



1 — 2(cotha; a ) 2 H cothx a 



(12) 



x a takes the same values as in Eq. (11). When x a — > D c (T,n a ) — * 1, how- 
ever this function decreases with x a much faster than D x (T,n a ). Therefore 
with decreasing uj a stronger damping of the oscillations in C is observed as 
compared with x- Besides, the amplitude of the dHvA oscillations in C is pro- 
portional to W 2 that additionally smoothes these oscillations with increasing 
e F /u. 

Another type of high frequency oscillations is connected with VF_-frequency 
through the function cos(2n \n^/W^). (Here frequencies must not be multiple 
again.) Thus the period of this oscillations rapidly falls down with increasing 
uj. Therefore the temperature damping becomes stronger at large uj. Hence one 
can find T start supposing that this temperature is such that D x (T start , nJ) ~ 0.1 
i.e. this type of oscillation is essentially suppressed at T star t. As D x = 0.1 at 
x ~ 4.5 and = Vt 2 /uj near the quantum limit one gains: 

T start ~ ^(3^)^. (13) 



In the 3D-case there could be exist one more type of oscillations: cos(27r-^-rio) 
(Wq/W + 7^ n /n+; W /W_ ^ n /n_). However these oscillations can really 

— * 

reveal themselves when the direction of B does not coincide with any of oscil- 
lator symmetry axes and only at small uj because near the quantum limit W 
does not depend on uj and these oscillations do not arise. 

In all three cases considered above we marked in brackets that W + , W_, W 
should not be multiple to each other. Nevertheless, if uj varies continuously 
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a pair of frequencies or even all three can become multiple. When u reaches 
such values that 

n + n_ n no n_ n + 

these types oscillations are suppressed at T = T start and uj < ep . The third type 
of such oscillations n + /W + = n /W ^ n^jW^ appears in such points of uj 
where the ratio of W + and W is equal to an integer, if n = 1, i.e. these points 
are divided by an interval approximately equal to W ~ Q. Therefore this type 
of oscillations can be damped by the temperature if _D x (T opt , n ) — 0.1. The 
parameter x in this case at n — 1 is 2n 2 T opt /W = 4.5: 

T opt ^i(3iV) 2/3 A. (14) 



The last type of oscillations W + /Wq/W- = n + /no/n_ is also damped because 
the minimum value of n_ is 1 and at uj > 2Q the frequency W_ <C Wq, i.e. 
riQ ^> 2 and D x <^ 0.1. Amplitudes of all oscillations with multiple frequencies 
increase with decreasing u and at small u < 2fl the temperature damping is 
neutralized. This is confirmed by our calculations indicating that in a wide 
range of N the cluster dHvA oscillations can be separated from others at 
T = T opt when the cyclotron frequency cj min is about two times larger than fl 
Eq. (5): 

ccw ^ 2SI, n max ~ ^iV) 1 / 3 . (15) 



Eq. (15) implies that these oscillations begin at such magnetic field, the cy- 
clotron radius R c of which turns out to be equal to the minimal size of a cluster 
in the plane perpendicular to the magnetic field. Thus u min is determined by 
the same conditions as in bulk metals. The existence of n max depending on iV 
limits the cluster particle number iV min at which the dHvA oscillations might 
be still observed. Since the minimum value of n max is 1 (the only parabolic 
segment shows up in x) Eq. (15) gives 7V min ~ 20 (Fig. 6). Fig. 6 and our cal- 
culations for iV ~ 10 4 show that for an array of clusters with wide enough size 
distribution ( and consequently with wide distribution in particle numbers) 
one can observe practically the same number of the oscillations as for a single 
cluster with N . 

We have analyzed the dependence of our results on variations of cluster shapes 
and the direction of the magnetic field. For this purpose at a fixed particle num- 
ber we have altered the cluster shape from oblate to prolate (0.1 < a z /a x < 3) 
with different degree of nonaxiality (1 < a y /a x < 2) and varied the direction of 
a magnetic field. Our calculations have shown that the properties of the dHvA 
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oscillations in 3D finite systems are rather insensitive to the variations of the 
cluster shape if the size of a cluster in the direction of the field is not too small 
as against the size in the perpendicular plane (e.g. a z should be greater than 
0.5(22,). Thus reasonably large variations of cluster shape and direction of the 
magnetic field cannot hinder measurements of the cluster dHvA oscillations. 

Values of minimal magnetic fields required to observe cluster dHvA oscillations 
at T opt can be assessed from Eq. (15) (obviously the upper limit is Sp) 

£f r (t , in4 m* e F (eV) 

^~2(aJvp' 5min(T) ~ 10 ^2^NW (16) 



These equations indicate that reasonable fields can be used for large grains 
of such materials for which the Fermi energy (sf) an d effective electron mass 
(to*) are rather small quantities. In addition, at a small Sp and large N the 
temperatures can be low enough, T opt ~ IK, to neglect electron scattering 
effects. For GaAs systems (m* = 0.067m, e F 10 meV) 

fi min (T) ~ 2/iV 1 / 3 ; T opt (K) ~ 20 /N 1 / 3 . (17) 



4 Summary 



In conclusion using the three-dimensional oscillator model we predict the 
dHvA oscillations in grains of materials with small Ef and to*. Analogously 
with bulk metals the cluster dHvA oscillations arise when the cyclotron radius 
is less than the size of a system. Size effects are clearly observed for N < 10 5 : 
the number of the oscillations is reduced with decreasing N; the size limit for 
N is ~ 20; the period of the oscillations stretches with increasing ef/uj. The 
properties of the dHvA oscillations (the number of oscillations, the amplitude, 
the period) are weakly sensitive to the significant variations of cluster shapes 
and the direction of the magnetic field especially near the quantum limit. The 
characteristic feature of the cluster dHvA oscillations is a special tempera- 
ture regime at which high frequency oscillations caused by the discreteness 
of electronic level spectra are suppressed. The oscillator model predicts that 
the 3-D-mesoscopic dHvA oscillations appear at T start ^ ef/3(3N) 2 ^ 3 and the 
most favorable temperature for the observation of their maximum number is 
T opt ~ e i ?/4(3A^) 1 / 3 , while an excess of T opt results in smoothing and subse- 
quent disappearing the cluster dHvA oscillations. 

We thank V.E. Bunakov for useful discussions. 
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Fig. 1. Low temperature oscillations of the electron susceptibility \ ( a ) an d 
heat capacity C (b) in a tilted magnetic field (uj z = uj cos 6, uj x = uj sin 8 cos(p, 
ujy = uj sin 6 sin ip, 8 = it/ 4, (p = vr/4) vs e_p/w for N = 10 5 electrons moving in an 
anysotropic oscillator potential (a y /a x = 1.33, a z /a x = 1.55). xl = — ^"bN/2Vsf, 
A 3.Y. 2. 
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Fig. 2. Oscillations of the magnetic susceptibility vs ef/uj for a cluster with N = 10 5 
electrons at T = T start = 22A. The oscillator anysotropy and the direction of the 
magnetic field are the same as in Fig. 1, g = 2. 
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Fig. 3. Influence of the temperature on the dHvA oscillations in a system with 
N = 10 5 . The oscillator anysotropy and the direction of the magnetic field are the 
same as in Fig. 1. 



12 




8 p / CO 

Fig. 4. Comparison of the dHvA oscillations in systems with N = 10 5 and N = 10 3 
at T = T opt . T opt = 10 3 A for N = 10 5 and T opt = 50 A for N = 10 3 . The oscillator 
anysotropy and the direction of the magnetic field are the same as in Fig. 1, g = 2. 
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Fig. 5. Influence of the effective Lande factor on the dHvA oscillations in systems 
with iV = 10 4 at T = T opt = 250A. The oscillator anysotropy and the direction of 
the magnetic field are the same as in Fig. 1. 
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Fig. 6. N-averaged dHvA oscillations of spherical clusters with N = Nq ± 3AN 
(Nq = 20, AiV is the averaging width parameter), g = 2. 
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